Dynamical Mean Field Theory for Perovskites 
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Using the Hubbard Hamiltonian for transition metal-3d and oxygen- 2p states with perovskite 
geometry, we present a dynamical mean field theory which becomes exact in the limit of large 
coordination numbers or equivalently large spatial dimensions D. The theory is based on a new 
description of these systems for large D using a selective treatment of different hopping processes 
which can not be generated by a unique scaling of the hopping element. The model is solved using 
a perturbational approach and an extended non-crossing approximation. We discuss the breakdown 
of the perturbation theory near half filling, the origin of the various 3d and 2p bands, the doping 
dependence of its spectral weight, and the evolution of quasi particles at the Fermi-level upon 
doping, leading to interesting insight into the dynamical character of the charge carriers near the 
metal insulator instability of transition metal oxide systems, three dimensional perovskites and other 
strongly correlated transition metal oxides. 



I. INTRODUCTION 

The electronic structure of the strongly correlated cop- 
per oxide superconductors, other transition metal based 
perovskites, and of transition metal oxides like NiO, 
FeO or MnO has been an enduring problem in the last 
years still containing numerous unresolved ques- 

tions. Among them are the nature of the insulating state, 
the evolution and doping dependence of coherent quasi 
particles near the Fermi energy in a doped Mott Hub- 
bard or charge transfer insulator Q and the transfer of 
spectral weight from high to low energy scales ||. Ex- 
perimentally, interesting variations of the spectroscopic 
and thermodynamic properties of these materials occur 
upon electron and hole doping P-|l2|]. Recently, Tokura 
et al. 1 10 1 showed that Sri-^La^TiOs with perovskite- 
like structure exhibits a pronounced increase of the ef- 
fective mass in approaching the insulating state while 
keeping the existence of a large Fermi surface which sat- 
isfies the Luttinger sum rule. They already argued that 
the different behavior of this compound compared to the 
pronounced two dimensional high-T c cuprates might be 
due its three dimensional character. Furthermore, the 
transfer of spectral weight in Li-doped NiO Jll]] and the 
investigation of hole and electron doping of this com- 
pound in Ref. fl^| demonstrate the importance of strong 
electronic correlations in transition metal oxides. A sys- 
tematic analysis of the electronic structure of 3d- transi- 
tion metal compounds was recently performed by Saitoh 
et al. |{|. They showed that for a variety of materials 
like LaMn0 3 , LaFe0 3 , LaCo0 3 , Mn .25TiS 2 , Ni . 33 TiS 2 , 
NiO, FeO, or MnO the character of the band gap is of 
the charge transfer type. Furthermore, the analysis of 
x-ray photoemission spectra [jb| demonstrates the im- 
portance of the charge transfer mechanism in systems 
like V2O5, T1O2 or LaCr0 3 and showed that many of 
the early transition metal compounds have to be reclas- 
sified as intermediate between the charge transfer and the 



Mott-Hubbard region or as falling into the charge trans- 
fer region itself. Consequently, it is of importance for an 
understanding of these materials that one takes the local 
spin and charge fluctuations of the transition metal and 
the oxygen states into account, i.e. both orbital degrees 
of freedom have to be described on the same footing. 

For systems with one orbital degree of freedom a theo- 
retical approach that allows the calculation of the excita- 
tion spectrum of strongly correlated systems, maintain- 
ing the dominating local correlations, was recently pro- 
posed by Metzner and Vollhard |Q. They introduced 
a dynamical mean field theory, where the system can 
be mapped on a local problem coupled to an effective 
bath. [p~5| — p~8|| This dynamical mean field approach is ex- 
act in the limit of large coordination numbers or equiv- 
alently for large spatial dimensions (D — ► 00). fl4| , |l9|| In 
most applications of this approach, the one band Hub- 
bard model is considered. A scaling of the hopping el- 
ement like t = t* 1 ' \f~D with fixed t* leads for D — > 00 
to a remarkable simplification of the many body prob- 
lem while retaining the nontrivial local dynamics of the 
elementary excitations and other main features of the 
model. 0-H 

Only few attempts have been made to extend the dy- 
namical mean field theory to systems with more than one 
orbital degree of freedom. Valenti and Gros |^4j realized 
that for the two band Hubbard model with perovskite 
structure, a scaling procedure similar to the one band 
model does not lead to a finite bandwidth. In order to 
avoid these difficulties Georges et al. introduced a 
CuO model where two interpenetrating hypercubic lat- 
tices of copper and oxygen sites have been considered. 
This leads to a finite bandwidth and is able to describe in- 
teresting physics of the model under consideration. How- 
ever, square root divergencies of the uncorrelated density 
of states occur, which result from the sublattice structure 
of the underlying lattice. This could limit the range of 
applicability of this model in the description of real tran- 
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sition metal oxides or possibly high-T c systems. 

In this paper we propose a scaling procedure to de- 
velop a reasonable model of perovskite systems for large 
D based on a selective treatment of different hopping pro- 
cesses in large dimensions. It will be shown that this pro- 
cedure leads to a large D version of the three-band Hub- 
bard model which reproduces important physics of the 
low dimensional situation. First, we solve this model us- 
ing perturbation theory (l7]]2(| and show that this leads 
to a qualitatively wrong description of the physics near 
half filling. These problems are avoided by using an ex- 
tended version of the non crossing approximation |27] [H| . 
Finally, we discuss the doping dependence of the transi- 
tion metal and oxygen densities of states. This is of inter- 
est due to the evolution of coherent quasiparticles near 
the Fermi energy, the transfer of spectral weight and the 
change of the effective mass upon doping. 



II. THEORY 

The dynamical mean field theory, which becomes exact 
in the limit of large coordination numbers, turned out to 
be a natural starting point for the description of highly 
correlated electronic systems like the one band Hubbard 
model and the Anderson lattice model. All this is based 
on a proper scaling of the hopping element, originally 
proposed in Ref . . In the following, we show that the 
extension of this approach to more complex systems is 
only possible if one generalizes the idea of the scaling pro- 
cedure in a way that the corresponding transfer functions 
between different unit cells, but not the hopping clement 
itself, have to be scaled. This is demonstrated for the case 
of the two band Hubbard Hamiltonian, which describes 
hybridized transition metal (TM) and oxygen states with 
perovskite lattice structure. Similar to the one band case, 
the problem can be mapped onto an Anderson impurity 
problem and the enormous amount of theoretical tools, 
developed for the solution of this model, can be used. 
Furthermore, we show that the most natural choice of 
the corresponding impurity problem is to include all lo- 
cal orbitals into the impurity, not solely the strongly cor- 
related ones. Otherwise, the set of self consistent equa- 
tions is not well defined, leading to various computational 
problems. In our case, a cluster of hybridized TM- and 
O-sites is embedded into an effective medium. Physically, 
this is necessary because the charge and spin fluctuations 
of these two states are closely intertwined and have to be 
described on the same footing. The solution of the ex- 
tended Anderson impurity model is performed using a 
perturbational approach [T^ , |2^ ] and an extended version 
of the non crossing approximation p7|-pT|. The details 
of the latter method are presented in the appendix. 



A. Scaling procedure for large D 

We consider the D-dimensional extension of the per- 
ovskite lattice, where the transition metal sites sit on a 
Z)-dimensional hypercubic lattice and the D oxygen sites 
per unit cell are located between every two nearest neigh- 
bor TM sites. For simplicity, we consider only one TM 
3d orbital. The corresponding (D + l)-band Hubbard 
Hamiltonian for the TM 3d and O 2p Xa (a = 1, . . . , D) 
orbitals reads 

1,(7 /,CT,a 

+ * E 9u(dtpia Xa +H.c.) + t/]T« ■ (1) 

il,(T,a i 

Here, d\ a (pl ax ) creates a hole in a TM-3d (0-2p Xa ) 
orbital at site i (I) with spin a. £d and e p are the cor- 
responding on-site energies, t is the amplitude of the 
nearest neighbor TM-0 hopping integral and U is the 
Coulomb repulsion between TM holes with occupation 
number operator nf a = d\ a di a . The p Xa — d hopping ma- 
trix elements have pdcr-symmetry, so that = +1 or — 1 
if I = i + x a /2 or I = i — x a /2, respectively. Only the 
combination p sia = -j= J2i a 9aPi<rx a of oxygen states 
hybridizes with the d-states p^j. Since the operators 
p s i a do not anticommute, it is necessary to orthogonalizc 
the corresponding states leading to pko- = 7k 1 Psk.a where 
Pskcr is the Fourier transform of p S i a and 

1 D 

7k = 1-— E cosfc «- ( 2 ) 

a=l 

Besides the Pko-, which now fulfill standard fermion com- 
mutation relations, D — 1 linear independent combina- 
tions of the Fourier transforms of p\ ax occur, which are 
orthogonal to p ko . and build up nonbonding oxygen or- 
bitals. These nonbonding states at e p decouple from the 
remainder of the system and do not have to be considered 
in the following calculation. The resulting Hamiltonian 
reads: 

+V2Dt J2 (ikdL-Pk" + H - c ) ( 3 ) 

ko- 

This is a straightforward extension of the two band Hub- 
bard model of perovskite systems to arbitrary dimension 
and can be used as the basis of a 1/D expansion. How- 
ever, there exists no scaling of the form t = t*/D* with 
exponent (3 leading to a nontrivial limit for D — > oo. 
Valenti and Gros [^4| suggested a scaling with = 1/2 
which leads to a zero bandwidth of the bonding and an- 
tibonding states, i.e. the TM-0 "dimers" of different 
lattice sites decouple. All other choices of (3 lead either 
to zero or to infinite density of states. 
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In order to obtain a physically reasonable limit of the 
perovskites for large D one has to bear in mind that dif- 
ferent hopping processes which are of the same order in t 
are of different importance for large D. Due to the local 
character of the Coulomb interaction, it is sufficient to 
perform the scaling at first for [7 = 0. For the calcu- 
lation of the single particle Green's function one has to 
sum up all closed paths on the lattice under considera- 
tion. In Fig. |l| we show two different hopping processes 
which are of forth order in t. There are AD 2 processes of 
type (a) and only 2D processes of type (b). Using a scal- 
ing like t — t*/y/D the contribution of the (a) processes 
to the Green's function are of order 1 whereas the (b) 
processes have only contributions of order 1/D and are 
suppressed for D — > oo leading to the zero bandwidth of 
Ref. f24|j . This results from the TM coordination num- 
ber being 2D whereas the oxygen coordination number is 
2. A selective consideration of these different paths can 
be obtained by the following partial summation of paths 
starting at TM-site i and ending at the TM-site j: 

G° dij (w) = GfWij +Gf{u)T ij {u)Gf{u) 

+ Y J Gt{u)T ll {u)Gf{u)T lj {u)Gf{u) +••• (4) 
i 

Here, the local d-Green's function 

Gf{u)=L-e d -^ r ) J (5) 

contains all TM-O-TM hopping processes, where the 
path returns to the starting TM-site without entering 
any other TM site ((a) processes). The transfer function 
Tij(w) which is given by (—t 2 )/(u> — e p ) for neighboring 
TM sites and zero otherwise has the character of an effec- 
tive, but frequency dependent TM-TM hopping element 
Tij ex 5(ij)t 2 ((b) processes). Therefore, we propose the 
following scaling procedure, where the two different hop- 
ping processes are treated differently: 

Gf(u) =L-e d - ; t* = ^D/2t , (6) 

T ij (a;) = -^ n /J-&; (f ) 2 = ^D/2t 2 . (7) 

V U LU — Ep 

Here, the effective TM-TM hopping (not t itself) has 
been scaled in analogy to the one band case. Conse- 
quently, the processes of Fig. |l|b are no more suppressed 
relative to that of Fig. [l]a. The summation of all closed 
paths on the remaining hypercubic TM-lattice is straight- 
forward and leads to the d-Green's function for U = 0: 

A{t*y J {uj - e p ){u - e d )/A(t*) 2 - y 
where the reduced density of states g {y) is defined by 



k 

with ak = X^=i cos ^- The oxygen Green's function 
G°(lu) which refers to the states Pk a can be obtained from 
G d (cu) by interchanging e d and e p . The large D behavior 
of g {y) is, in analogy to the one band case, given by 
a Gaussian distribution function ^/2/7rexp(— 2y 2 ). |ll]] 
For simplicity, we use in the following a semi elliptical 
reduced density of states g (y) = 2/ir^l — y 2 if y 2 < 1 
and zero otherwise. This Bethe lattice of the TM-sites 
does not changes the relative local arrangement of TM 
and O sites, typical for perovskite systems. Finally, it 
is interesting to point out, that this selected treatment 
of different hopping processes can be generated by the 
Hamiltonian of Eq. if one performs the following scal- 
ing transformation of the coherence factor: 

t 2 D 1 2 ^2(t*) 2 ~2^(n 2 a k . (10) 

Based on this formulation of the large D version of the 
TM-0 system for U = 0, it is straightforward to show for 
U 7^ Q, in analogy to the one band case, that the self en- 
ergy Sd(w) of TM-holes is momentum independent. Con- 
sequently, one finds for the TM-Green's function within 
the Bethe lattice: 

4i 2 

G d (uS) =lj- e d - 

LU - Ep 

1 / At 2 \ 2 

4 \lj -s p J 

This results from the representation of the Green's func- 
tion with semi-elliptical density of states: 

k 



uj -e p ( (cj - E p ){ui -e d - Srf(^)) 
At 2 \ At 2 




with T{z) = 2 yz — 1 — \J (z — l) 2 — 1 j and the relation 

J-(z)^ 1 = z — 1 — jJ 7 (z). Furthermore, the system can 
be mapped onto an Anderson model with effective hy- 
bridization and additional self consistency condition. 

Finally, we note that the additional consideration of 
the TM-0 Coulomb repulsion U pd would lead to fur- 
ther frequency dependent contributions of the self ener- 
gies E Q| a(w) (a/3 € {p, d}). This results from the small 
coordination number of the oxygen sites and is different 
in the one band case, where each intersite Coulomb re- 
pulsion contributes in zeroth order of 1/D only through 
its Hartree-Fock value Mj. Therefore, the limit D — > oo 
might be an interesting approach for a detailed consider- 
ation of the dynamical interplay of spin and charge fluc- 
tuations of transition metal oxides. 



3 



B. Effective TM-O impurity model 

As in the one band case, the large-D version of a TM-0 
model is mapped onto an effective impurity model. Due 
to the occurrence of an additional orbital degree of free- 
dom, the choice of the corresponding effective medium 
is not unique. Although, the most obvious choice might 
be to embed only the TM-site into an effective medium 
and to change solely the self consistence condition, it is 
straightforward to show that this leads to a divergence 
of the hybridization of this medium at the oxygen on- 
site energy e p . This is due to the fact that the effective 
medium has to simulate not solely the correlated TM- 
sites but also the existence of the oxygen states. In order 
to avoid this divergency, we propose in the following an 
impurity model where a cluster of one TM- and one O- 
orbital is embedded within an effective medium. This is 
illustrated in Fig. ^. Motivated by the perovskite geom- 
etry, we couple the effective medium only to the oxygen 
sites. Furthermore, the local TM-0 hybridization within 
the cluster is explicitly taken into account. This leads 
to a well defined effective medium. The corresponding 
Hamiltonian reads: 

H e s = H\ oc + H mc d , (13) 

where the local part is given by 



and 



kcr 



ktr 



describes the coupling of the oxygen states with the ef- 
fective medium, characterized by the hybridization 



U3 + l0 + — £k 



(14) 



This coupling of the effective medium to the oxygen 
states, although not unique, is physically motivated by 
the local arrangement of the original lattice. From the 
equation of motion of the effective Hamiltonian, it follows 



G d (u)) 1 = w - e d - 



Comparing this with Eq. [ll], one immediately finds 

t 2 G d (uj) 
uj - £ p - t 2 Gd{w) 



j(u) = rG d (cj) i - 



(15) 



(16) 



This is the self-consistant equation of the theory, based 
on the condition that the Green's function of the impu- 



rity model defined in Eq. 13 equals that of the lattice 



system. Note that this effective medium, which repre- 
sents the states coupled to the oxygen, is predominantly 
determined by the correlated TM-sites. 

In the following we solve the impurity model of Eq. [l3| 
using the iterated perturbation theory (IPT), and an ex- 
tended version of the non crossing approximation (NCA) , 
where the effective medium is fixed by Eq. [l6|. Since it 
was necessary for the consideration of the TM-0 clus- 
ter to generalize the standard NCA-scheme to more local 
eigenstates, the corresponding details are given in the 
appendix. 



III. RESULTS 

In the following, we present our results for the TM- 
and O-densities of states obtained within the IPT and 
NCA. In Fig. [| we show our results for the TM- and O- 
densities of states obtained within the IPT for half filling 
x = (x — n d + n p — 1) . Two Hubbards bands which 
are dominated by TM-states are separated by U and two 
oxygen dominated bands in the neighborhood of e p are 
visible. Similar to the one band case, the IPT leads to 
an interesting structure of the density of states, where 
now four bands are occurring. However, the expected 
insulating behavior at half filling (for large U) does not 
occur. This results from the overestimation of the spec- 
tral weight of the lower Hubbard band within the IPT by 
cx 1 — rid- The occupation of oxygen sites (n p > 0) due to 
TM-0 hybridization leads to an overcounting of copper 
sites which can be occupied without paying any Coulomb 
energy. Therefore, the success of the IPT in the one band 
case |l7j] does not occur within the two band model. This 
is due to the absence of the particle-hole symmetry and 
the change of the TM-occupation number rid as function 
of t at half filling. This shortcoming of the IPT results 
from the occupation of O sites and is expected to vanish 
for small n p which occurs in the limit of a large value for 
A = e p — Ed- However, for the physically interesting situ- 
ation A k, U/2, the IPT leads to qualitatively wrong re- 
sults and one has to develop theoretical approaches which 
take the local TM-0 many body states explicitly into ac- 
count. Therefore, we use the NCA , which was shown 
to be in excellent agreement with quantum Monte Carlo 
(QMC) simulations in the one band case, such that a lo- 
cal TM-0 hybrid with 16 local eigenstates is coupled to 
an effective medium. Here, the TM-0 singlet state which 
is related to an important excitation near half filling is 
explicitly taken into account. 

In Fig. ^ we show the TM- and O-densities of states 
at half filling and for a hole and an electron doped sys- 
tem, obtained within the non crossing approximation. 
The calculations are performed for U = 2A = lOt and 
t = 1 eV. At half filling we find now the expected insu- 
lating state, i.e. the Fermi level lies in the middle of the 
charge transfer gap. In this charge transfer regime we 
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find bands dominated by TM-states (B, F, and G) which 
are separated by the Coulomb repulsion U due to the 
Mott-Hubbard splitting and on both sides of the charge 
transfer gap TM- and O-dominated bands (B and C). 
Besides the O-dominated band above the charge transfer 
gap (C-D) a second one for higher excitation energy oc- 
curs (E). A detailed analysis of the transitions between 
the various many body states of the local TM-0 hy- 
brid reveals that these are build up by the TM-0 sin- 
glet and triplet states. All these features of the den- 
sity of states are believed to be general properties of 
the model under consideration [ [33] and are also quali- 
tatively in agreement with exact cluster diagonalizations 
and Monte Carlo simulations for two dimensional sys- 
tems. [^i|-^6| This demonstrates that it is indeed possi- 
ble to reproduce important physical phenomena of the 
low dimensional situation within the large D approach. 
Besides this general structure on a larger energy scale, a 
coherent quasi particle peak near the Fermi level occurs, 
which shows up an interesting temperature and doping 
dependency (see below). 

In order to clarify the origin of the details of the den- 
sity of states, we investigated which transitions between 
the local many body states are predominately responsible 
for each band. Technically, this is performed within the 
NCA if one restricts the available Hilbert space according 
to the general scheme presented in the appendix. If the 
spectral weight of a band decreases after projecting out a 
certain local state |m), a transition from or into this state 
is responsible for the band. This analysis is of particu- 
lar importance, if one is interested in calculations within 
a restricted Hilbert space of the most important states. 
The results are presented in Fig. ^) for hole doped sys- 
tems, whereby the meaning of the state labels are given 
in Tab. |. The leading contribution is underlined while 
negligible contributions are embraced. 

The main contributions to the bands B and E occur 
also for an uncorrelated system and are build up by tran- 
sitions between an empty state |0) and a bonding (\d a )) 
or antibonding combination of TM- and O-states, 

respectively. More interesting are the bands A, C, F, 
and G, because they are purely due to strong electronic 
correlations. Here, the bands F and G are due to transi- 
tions from the TM-dominated state \d a ) to states which 
are dominated by doubly occupied TM sites, i.e. they 
form the upper Hubbard band. Band C contains essen- 
tially the well known TM-0 Zhang Rice singlet state J3^] 
resulting from a transition from the bonding state \d a ) 
to the state \S Q ) dominated by the TM-0 singlet state, 
other contributions close in energy are identified under la- 
bel D. Furthermore, a pronounced part of band E, which 
consists predominantly of antibonding \p a ) states, results 
also from the transition of singly occupied sites to TM- 
O triplet states \T±) . Consequently, we find a singlet 
triplet splitting which compares well with the result 8J 
obtained from perturbation theory in t for two spatial di- 



mensions, where J = t 2 {hr l + (U - A)" 1 ). |37) Finally, 
the lowest state A results from a transition between the 
antibonding state \p a ) and the singlet state \S ). Inter- 
estingly, this state occurs for very low excitation energies 
and only in systems with an occupied singlet state, i.e. 
for hole doped systems. This is a consequence of the fact 
that the singlet state with two particles per TM-0 unit 
has lower energy than the singly occupied antibonding 
state \p a }- 

This interpretation of the various bands in terms of the 
transitions between many body states allows a transpar- 
ent explanation of the doping dependence of their spec- 
tral weight. In Fig. || we show the spectral weight of 
various bands upon electron and hole doping. For hole 
doping we plot in part a of the figure the spectral weight 
of all occupied states and of the occupied part of the 
singlet band (low energy hole states). The latter demon- 
strates that the transfer of spectral weight, discussed in 
Rcf. is included within our approach. As expected, 
the spectral weight of the bonding band decreases upon 
doping. However, band A which results from the transi- 
tion between \S ) and \p G ) increases, because upon hole 
doping more initial singlet states are available The cor- 
responding result for the electron doped materials are 
shown in Fig. Ifa, where we plot the spectral weight of 
all unoccupied bands and that of the unoccupied part of 
the lower Hubbard band B which carries the Fermi level 
(low energy electron states) . While the spectral weight of 
the bands which are purely due to correlations decreases, 
band E increases upon electron doping. The character of 
the system becomes much more uncorrelated leading to 
an increase of the bonding and antibonding TM-0 states 
which are present also for U — with even larger spectral 
weight. 

Based on this analysis, which demonstrates that the 
density of states near the Fermi level is build up by the 
transitions between only a few local many body states, 
it is interesting to perform the calculation within a re- 
stricted Hilbert space. In Fig. || we show the TM- and 
O-densities of states if one takes only the states |0), 
\p a ) ,\d a ) ,\S ) , and \T±) into account. The results are 
compared with that for the full Hilbert space (dashed 
line). It is impressive, that the details of the single par- 
ticle excitation spectrum near the Fermi level are very 
well reproduced by this restricted set of states. Note that 
this calculation includes still the full information about 
the orbital character of the low energy spectrum. This 
might be of interest for the theoretical description of more 
complex systems where the diagonalization of the local 
problem leads to extremely large numbers of coupled in- 
tegral equations of the NCA. The most important states 
can straightforwardly be found from the eigenstates of 
the complete local problem, since the peak position is 
well reproduced from the corresponding difference of the 
relative eigenenergies, as given in Tab. [j| 

Finally, we present our results for the coherent low 
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energy peak near the Fermi level. In Fig. [7] we show 
the density of states near the Fermi level (u> — 0) for a 
hole doping concentration x = 0.2 for various temper- 
atures |3§| ]. For electron doping, the corresponding low 
energy part of the spectrum behaves very similarly. Al- 
though the density of states at the Fermi level is almost 
constant, the states near the Fermi energy are strongly 
temperature dependent. For decreasing temperature, a 
coherent quasi particle peak occurs which builds up its 
maximum slightly below the Fermi energy. 

In Fig. ||a and b we present our results for the dop- 
ing dependence of the low energy part of the spectrum 
for a hole doped system in the charge transfer and for 
an electron doped system in the Mott-Hubbard regime 
at T = 300 K. As expected, the general behavior of the 
low energy part of the total density of states is relatively 
similar in both regimes. The evolution of these quasipar- 
ticles results in a large effective mass ratio obtained via 
m* /m — (1 — 9ReSrf(a;)/9a;)| w=0 . In Fig. ||c, as a rep- 
resentee example we show results for a electron doped 
material in the charge transfer regime. This shows clearly 
the mass enhancement when approaching the insulating 
regime either charge transfer or Mott Hubbard type and 
resembles the doping dependence of the effective mass 
obtained within the Brinkmann-Rice transition 
where only the coherent part of the spectrum occurs. In 
addition, we plot in part a of Fig. || the density of states 
for low hole doping (x = 0.017) but for large tempera- 
tures of T = 1000 K. This shows that the quasi particle 
peak disappears completely between 300 and 1000 K, giv- 
ing rise to an enormous temperature dependence of the 
low energy excitations in slightly doped perovskite sys- 
tems. The enhancement of m* /m with band filling agrees 
with the observed behavior in Sri-^La^TiOa [[10| which 
is believed to fall into the the Mott-Hubbard regime. Fur- 
thermore, the measurement of the Hall coefficient indi- 
cates the formation of a large Fermi surface. Although 
our dynamical mean field approach does not allow an in- 
vestigation of the momentum dependent spectral density 

A(k, u) « —6{lu - 4) + A incoh - (k, w) , 
m 

we believe that a large effective mass ratio m* jm leads 
to a large Fermi surface formed by the coherent part of 
the spectrum ^8(w — e£) with renormalized quasi par- 
ticle energy e£. This is possible because only a few par- 
ticles can occupy a large number of momentum states 
due to its low spectral weight j^-. The physical origin 
of the coherent quasiparticle resonance is similar to the 
conventional Abrikosov-Suhl resonance of systems with 
magnetic impurities. In the charge transfer regime, for 
hole doping, the resonant transition between a bonding 
state \d a ) and the Zhang Rice singlet \S ) builds up the 
Kondo like peak. For electron doping, it is correspond- 
ingly a transition between an empty state |0) and \d a ). 



IV. CONCLUSIONS 



In conclusion, we presented a dynamical mean field 
theory for perovskites which is based on a new scaling 
procedure of the large dimensional Hubbard model of 
transition metal 3d and oxygen 2p states. In the limit 
D — > oo, the electronic self energy is, similar to the one 
band case, momentum independent and the model can be 
mapped onto an Anderson impurity model. It is demon- 
strated that the most natural choice of this impurity is 
to include all local orbitals of one unit cell, not solely 
the strongly correlated ones. Correspondingly a model 
with 16 local many body states is embedded into an ef- 
fective medium which has to be determined self consis- 
tently. Motivated by the perovskite lattice geometry, this 
effective medium is coupled only to the oxygen sites of 
the impurity and simulates predominantly the TM-sites 
of the real lattice. The solution of the local problem is 
performed using the iterated perturbation theory and an 
extended noncrossing approximation. Is is shown that 
the spectral weight of the various bands of the density 
of states is erroneously described within the IPT. This 
shortcoming results from the missing electron hole sym- 
metry of the half filled model, leading to an absence of a 
metal insulator transition for half filling. In distinction, 
the NCA is able to describe this problem properly and 
gives interesting insight into the origin and doping depen- 
dence of the various bands of the TM- and O-densities of 
states. Furthermore, the occurrence and doping as well 
as temperature dependencies of a coherent quasi particle 
peak near the Fermi level is discussed. This results from 
a Kondo like resonance between an empty and a bonding 
TM-0 state for electron doping and this bonding state 
and the Zhang Rice singlet state for hole doping. Finally, 
it is shown that the restriction of the local states to a lim- 
ited set of most important states reproduces very well the 
low energy density of states, keeping the full orbital in- 
formation of the involved states. This is of importance 
for an extension of the theory to more realistic situations 
including all five d-bands such that only a small number 
of coupled NCA equations have to be solved. This might 
be stimulating for further theoretical studies of transi- 
tion metal compounds within the dynamical mean field 
approach. 
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APPENDIX A: NON CROSSING 
APPROXIMATION OF A TM-O CLUSTER 

In this appendix, we discuss an extended version of the 
non crossing approximation for an impurity model given 
in Eq. [t3| The NCA was shown to be in excellent agree- 
ment with quantum Monte-Carlo calculations for the dy- 
namical mean- field theory of the one band case p2| , p3|] . 
This approximation is not a perturbational expansion 
with respect to the Coulomb repulsion U, but takes into 
account all eigenstates of the local part of the impurity 
Hamiltonian, by performing an expansion with respect 
to the hybridization J(u>). Therefore, we rewrite the 
Hamiltonian of Eq. [l3| in terms of the Hubbard operators 
X m i m = \m')(m\, where \m) is, in our case, one of the 16 
eigenstates of iJ; oc , with eigenvalue E m (see Tab. [jj): 

P^loc — ^ -^m A mm ; 
am 

H m ed = ( W * C i<T P mm' X rnm> + H.C^j + ^ ^ c la C ^ 

Icamm' kcr 

The eigenvalues of Hi oc are explicitly given in Ref. p2| . 
The (16 x 16) matrix P a and a corresponding matrix D a 
are defined by the following representation of the O- and 
TM-destruction operators: 

Pa = P-mm' Xmm' , 

m,m' 

da — ^ ] Dmm'Xmm' ■ 
m,m' 

Using this representation of the Hamiltonian, one intro- 
duces within NCA scheme the local propagators for each 
of the eigenstates |m): 



P m (u) = 



1 



lo + i0+ - E m - £ m (w) 



(Al) 



with the corresponding spectral density p rn (uj) — 
— — ImP m (w + i0 + ). The NCA self energies are given 
by- 



1 f 

E m (w) = Y \P m ,m'? / def(r] m . m >£) 

lmj(e + iO + )P m > (w + rjm,m>e) ■ (A2) 

Here, f{e) is the Fermi distribution function and r\ m ,m' — 
+1 (rcsp. —1) if the particle number of the state \m') 
is higher (resp. lower) than the particle number of the 
state \m). In Eq. |A2|, we neglect the vertex corrections 
discussed in Refs. |31|]39| -|4l|]. Although this was shown 
to slightly underestimate the Kondo temperature and to 
shift the Abrikosov-Suhl resonance slightly away from the 
chemical potential, it gives important insights into the 
low energy excitations, in particular for large values of 



the Coulomb repulsion. Furthermore, in view of the enor- 
mous numerical problems, which arizes if one calculates 
these vertex corrections even for less complex systems 
than ours, [ fl0|]4l| ] and due to the additional self consis- 
tent condition of the large D approach, it is actually very 



difficult to go beyond Eq. A2. Finally, the Green's func- 



tions of the system, which give the information about the 
single particle excitation spectrum, result from: 

G d (u) - Yl \D mm '?((X m>m ;X mm ,)) u , (A3) 

m,m f 

G p (uj) = ]T \P^ n ,\ 2 ((X m , m ; X mm ,)) u , (A4) 



where ({X m > m ; X mm i)) u} is the Fourier transform of the 
retarded Green's function — i6(t){[X m r m (t); X mm r(Q)] + ). 
They can be expressed with respect to the NCA local 
propagators as following: 

1 f°° 

({X m > m ;X mm >)) u = — / dee~ Pe [p m -{e)P m {uj + e) 
-Pm(e)iV(w-e)*) , (A5) 

where the partition sum Zo is defined by 
Em /- dee~ l3e p m (e). All the summations with respect 
to m and ml are performed over the 16 eigenstates of 
Hi oc . As can be seen from Tab. § the spin degeneracy 
leads to only 11 different eigenstates and propagators 
in the paramagnetic state. The solution of Eq. Al and 



A2, which are coupled singular integral equations, is 
performed using the fast Fourier transformation for the 
convolution integrals and using the defect propagators of 
the local eigenstates [E8Le9|. 
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state label 


state 


degeneracy 


occupation n 


_E m - fin 


1 


10} 


1 








2 


\d a ) = (cos <j>dl + sin <j>pl)\0) 


2 


1 


-0.701 


3 


p CT ) = (cos^d* - sin^pt ) 0} 


2 


1 


5.701 


4 




1 


2 


5 


5 


|T CT )=dtpt]0) 


2 


2 


5 


6 


\S ) = [uidjdj + f 


1 


2 


2.738 




-«ji(4pI-dI P J)]|o) 








7 


\dd) = [uvdtd] + V2p\p\ 


1 


2 


12.217 




-w 2 (d\p\ -d\p\)]\Q) 

j. 4- 1 ' 








8 


\pp) = [u 3 dld\ + v 3 p\p] 

\r r / L [ 4- 14- 


1 


2 


10 




-w 3 (d\p\ -d[p\)]\0) 








9 


\D a ) = (cos 6dipU- a 










+ sin6dl P y_ a )\Q) 


2 


3 


9.298 


10 


\P a ) = {cobBcCIpIpU 










- sin 9dipUl<,)\a) 


2 


3 


15.702 


11 


\DP) =d\p\d\p\\0) 


1 


4 


20 



TABLE I. The degeneracy, occupation and relative energy of all local eigenstates introduced in the appendix and the 
corresponding labels used in Fig. 4 are given. The eigenvalues are obtained for U = 10 eV, A = 5eV, and t — leV. Here, 
tan 20 = — and tan 26 = -^jj ■ The coefficients m, Vi, and Wi can be obtained from the diagonalization of the three 
dimensional space build up by \S ), \dd), and \pp), see for example Ref. ^| . 
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FIG. 1. Hopping processes of order t 4 which are of different importance for large D. Since there are 4D 2 processes of type 
(a) but only 2D processes of type (b), the latter are suppressed compared to the former for large dimensions. 



N J(ffl) J(W)/ 




FIG. 2. Schematic impurity model where a cluster of one TM- and O-site is embedded within an effective medium. Here, 
the local hybridization t between the two local sites is explicitly taken into account. The effective medium is coupled only to 
the oxygen states. 
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FIG. 3. TM (solid line) and O (dashed line) density of states obtained within the iterated perturbation theory for half filling. 
Note that the spectral weight of the various bands is erroneously described within the IPT such that no insulator occurs for 
half filling. 
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FIG. 4. TM (solid line) and O (dashed line) density of states obtained within the hybrid-NCA at half filling (a) for hole 
doping x = 0.1 (b) and electron doping x = —0.1 (c). For hole doping it is indicated which transitions between the local states 
contribute to the corresponding bands, where the meaning of the state labels is given in Tab. 1. The leading contribution is 
underlined while negligible contributions are embraced. 
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FIG. 5. Doping dependence of the spectral weight multiplied with the spin degeneration of the occupied bands and the 
occupied part of the singlet band (low energy hole states) for hole doping (a) and of the unoccupied bands and the unoccupied 
part of the lower Hubbard band (low energy electron states) for electron doping (b) . 
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FIG. 6. TM (a) and O (b) density of states obtained within the hybrid-NCA for hole doping x = 0.1 within a restricted 
Hilbert space of the local states |0), |d CT ), \p a ), \T a ), and \S ) in comparison with the results of the full local space (dashed 
lines). 



0.4 



0.3 

^ 0.2 
? 

0.1 



0.0 





0)[eV] 











T=300K 








T=500K 








T=600K 








T=800K 








T=1000K 




r 1 


h 

%\ 



FIG. 7. Density of states near the Fermi level for hole doping x = 0.2 and for various temperatures. Note the formation of 
a coherent quasi particle peak resulting from Kondo like resonances between the singly occupied bonding TM-O state and the 
doubly occupied singlet state. 
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FIG. 8. Doping dependence of the coherent quasi particle peak for T = 300 K in the charge transfer regime with U = 10 eV, 
A = 5eV (a) and in the Mott-Hubbard regime with U = 5eV, A = 9.5 eV (b). In both cases we use t = 1 eV. The effective 
mass in the charge transfer regime and for electron doping is shown in part (c). For comparison we also plot in the inset 
the result for low hole doping (x = 0.017) but high temperature T = 1000 K which demonstrates that strong temperature 
dependence occurs in the low doping regime. 
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